Analytical Studies on a Modified Nagel-Schreckenberg Model with the Fukui-Ishibashi 

Acceleration Rule 

Chuan-Ji Pu, Chuan-Yang Yin, Tao Zhou, Bo Hu, Bing-Hong WangQ and Kun Gao 

Nonlinear Science Center and Department of Modem Physics, 
' University of Science and Technology of China, Hefei, 230026, PR China 

Q ! (Dated: February 2, 2008) 

. We propose and study a one-dimensional traffic flow cellular automaton model of high-speed ve- 

hicles with the Fukui-Ishibashi-type (FI) acceleration rule for all cars, and the Nagel-Schreckenberg- 
type (NS) stochastic delay mechanism. By using the car-oriented mean field theory, we obtain 
a : analytically the fundamental diagrams of the average speed and vehicle flux depending on the ve- 

hicle density and stochastic delay probability. Our theoretical results, which may contribute to the 
exact analytical theory of the NS model, are in excellent agreement with numerical simulations. 
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PACS numbers: 64.60.Ak, 89.40.+k. 02.60.Cb 



I. INTRODUCTION 



Traffic flow that displays various complex behaviors is a kind of many-body systems of strongly interacting vehicles. 
One of the approaches to microscopic traffic processes is based on cellular automaton (CA) Q], which treats the 
motions of cars as hopping processes on one-dimensional lattices. In the past few decades, CA models for traffic 
flow have attracted much interest of a community of physicists 0,0- Compared with other dynamical approaches, 
for instance the fluid dynamical approach, CA models are conceptually simpler, and can be easily implemented on 
'■^ | computers for numerical investigations 0. 

Two popular one-dimensional (ID) CA models are the Nagel-Schreckenberg (NS) model dj and the Fukui-Ishibashi 
(FI) model [5j, where periodic boundary conditions are used to mimic the traffic flow on highway. An exact car- 
, oriented mean field theory (COMF) have been developed to study the FI model, with an arbitrary limit on the 

maximum speed v max , car density p, and the delay probability / [fj. So far as we know, however, there has been no 
CsJ , established exact analytical theory for the NS model, due to the complications in the time evolution of the flow caused 
by the acceleration and stochastic delay rules. In order to comprehend how these rules affect the evolution and the 
corresponding asymptotic steady state, we study a ID traffic flow CA model which combines the NS stochastic delay 
£Tt ' rule and the FI acceleration mechanism. 

| This paper is organized as follows. In Sec. II, we give the definition of the model and the evolution equations for 
the inter-car spacings. Next, we present the fundamental diagrams of the average speed and vehicle flux depending 
' on the vehicle density and stochastic delay probability, compared with the numerical simulations. In Sec. IV, we 
f^) , summary with a discussion of our results in connection with the FI and NS models. 



II. THE MODEL AND ANALYTICAL SOLUTION OF ASYMPTOTIC AVERAGE VELOCITY 



The modified Nagel-Schreckenberg model with the Fukui-Ishibashi acceleration rule is a probabilistic automaton, 
q in which space and time are discrete and hence also the velocities. The road of length L is divided into cells of certain 
O ■ length, each of which can either be empty or occupied by just one car. The state of the nth car (n = 1, . . . , N) is 
characterized by the momentary velocity v n (t); in this model, the maximum speed is fixed as v max = M, and the 
vehicle density is p = N/L. Let C n (t) represent the number of empty sites in front of the nth vehicle at the time step 
t, then we have 

C n (t + 1) = C n (t) + v n+1 (t) - v n {t). (1) 

As a function of the inter-car spacing C n (t) and the stochastic delay probability /, the velocity of the nth car at time 
step t can be written as: 

v n (t) = F M (f,C n (t)), (2) 
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where 

0, if C = 

C - 1 with probability /, if < C < M 

C with probability 1 - /, if < C < M (3) 

M - 1 with probability /, if C > M 

M with probability 1 - /, if C > M. 

Here, as one can see, we have approximately adopted the acceleration rule introduced by Fukui and Ishibashi. The 
FI model can be considered as a NS model for "aggressive driving" , since the rules of the FI are nearly identical to 
the NS model, except that the acceleration rule has been changed from "the vehicle speed is at most increased by 
1 at each step" to "every car accelerates to v max " (which is adopted in the present model) and the stochastic delay 
mechanism has been modified as "only cars with v max will delay stochastically" (which is not adopted here). For 
v m ax = l, this may not change anything, however, for higher velocities it will lead to a considerable enhancement of 
the flow. Although the model is less realistic than the NS, it is still of interest due to its simplicity. The analytical 
description of the present model is of course much simpler than that of the NS, and hence it might serve as a testing 
ground for new analytical methods as well as new models. 

If Ni(t) represents the number of inter-car spacings with length i at time t, then the probability of finding such a 
spacing at time t is Pi(t) = Ni(t)/N. Hereafter, Pi(t) will indicated by P, for simplicity, except if specified otherwise. 
Let Qi denote the probability that an arbitrary vehicle moves i sites during a given time step, then we have: 

Qo = Po + fPi 
Q l = Pj(l - /) + /P i+ i(l <i< M -2) 

Qm-i = Pm-i(1-/) + /E£=m^ () 
Qm = (1-/)££"mP*. 

Obviously, Pi — when i > L — N, since the maximum inter-car spacing must be less than L — N. To obtain the 
nonvanishing Pi , we introduce Wi^j to denote the probability of finding an inter-car spacing with length i at a given 
time which changes into length j at the next time. According to the Eqs. (l)-(4), one can write all the nonzero Wi^j 
as follows: 



F M {f,C) = I 



W ^ j = P Q j , (l<j<M) 
W^o=m-f)Qo, (1<*<M-1) 
Wi^j = PiifQj-i + (1 - f)Qj], (1 < i < M - 1, 1 < j < M, j) 

W^m+i = PifQ M , (1 < i < M - 1) . . 

W^-m = P(l - f)Qo, (*>M) [0> 
W^i^M-j) = PlfQj-i + (1 - f)Qj], (i > M, 1 < j < M - 1) 
W^+i - PifQ M , (i > M) 
Wi^j = 0, otherwise. 

When the system approaches its asymptotic steady state, all the Pj will cease to change; thus the detailed statistical 
equilibrium condition for the steady state holds: 

]T = ]T VF m _,, Vm. (6) 

Since there are L — N + 1 variables Po, Pi, • • • Pl-n in Eqs. (6), which satisfy the following two normalization 
equations: 

L-N L-N 

Y, p k = i, E kPk = - - 1 ( ? ) 

fc=0 fe=l p 

we should truncate Eqs. (6) and obtain other L — N — 1 equations. In the present paper, we will choose the L — N — 1 
equations with m from to L — N — 2. Thus, we will have L — N + 1 independent equations, from which one can 
readily work out the values of Po, Pi, • • • Pl-n, then the average speed of the traffic in the steady state is: 

(v(t oo)) = Efli Pi[(i 1)/ + i(l - /)] + Ef=M+i Pi[(M - 1)/ + M(l - /)] 

= Efli iPi + Ef=M+i MP - f(l - Po) K> 
Based on the above discussions, one can easily obtain the traffic flux of the steady state: 

J(t-> oo) = p(v(t ->■ oo)). (9) 
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FIG. 1: The fundamental diagram of the average speed with the maximum speed M — 2 and for different stochastic delay 
probabilities /. The solid curves are numerical simulations while the points with different symbols represent theoretical results. 
The curves from the top down along the vehicle velocity axis correspond to different values of / ranging from to 0.9, in steps 
of O.f. 




FIG. 2: The fundamental diagram of the vehicle flux with the maximum speed M = 2 and for different stochastic delay 
probabilities /. The solid curves are numerical simulations while the points with different symbols represent theoretical results. 
The curves from the top down along the traffic flux axis correspond to different values of / ranging from to 0.9, in steps of 
O.f. 



III. NUMERICAL SIMULATIONS IN COMPARISON WITH THEORETICAL RESULTS 

In order to compare with the analytical results, we performed numerical simulations on a ID CA with length 
L = 1000 and the maximum car velocity M = 2. The number of vehicles N is adjusted so as to give the desired 
vehicle density p. The first 5000 time steps are excluded from the averaging procedure so as to remove the transient 
behavior. The average values are taken over the next 1000 time steps. Figure 1 and 2 show comparisons between 
results obtained from numerical simulations and car-oriented mean-field theory over the entire range of the vehicle 
density p. The curves are the simulation results, while the symbols represent the theoretical results. Theoretical 
results are in excellent agreement with numerical simulations. 
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IV. DISCUSSION 

In summary, we propose and study a one-dimensional CA model of high-speed vehicles with the FI acceleration rule 
and the NS stochastic delay mechanism. The analysis of the dynamical evolution of our model give us a clearly physical 
picture of how the acceleration and stochastic delay rules affect the evolution and the corresponding asymptotic steady 
state. Although cellular automata are designed for efficient computer simulation studies, an analytical description is 
possible, though difficult all too often. The method presented here can also be applied to other CA models. Our study 
may shed some new light on developing analytical approaches to other one-dimension traffic flow models, especially 
the NS model, for which until now, no exact analytical approach has been established. The investigation of the NS 
model has lead to a better understanding of its advantages and limitations, so that it is easier to choose the approach 
most suitable for a given problem. The question, how the stationary state is approached, is still an important open 
issue, which is currently under investigation and promises to yield new insights into the physics of the NS model. 
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